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Abstract. In a recent paper [H| J. Haglund showed that the expression (A^ . E n ^, e n ) 
with A hj the Macdonald eigen-operator A^if M = h^B^H^ enumerates by t area q dmv the 
parking functions whose diagonal word is in the shuffle 12 - ■ ■ jUJj + 1 ■ ■ ■ j + n with k of 
the cars j + 1, . . . , j -+- n in the main diagonal including car j + n in the cell (I, I). In view 
of some recent conjectures of Haglund-Morse-Zabrocki [12] it is natural to conjecture 
that replacing by the modified Hall-Littlewood funtions C P1 C P2 ■ ■ ■ C Pk 1 would 

yield a polynomial that enumerates the same collection of parking functions but now 
restricted by the requirement that the Dyck path supporting cars j + 1, . . . , j + n hits the 
diagonal according to the composition p = (pi,p2, . . . ,Pk)- We prove here this conjecture 
by deriving a recursion for the polynomial (Ah j C pl C P2 ■ ■ ■ C Pk 1 , e n ) then using this 
recursion to construct a new dinv statistic we will denote ndinv and show that this 
polynomial enumerates the latter parking functions by t area q ndlnv 



1. Introduction 



Parking functions are endowed by a colorful history and jargon (see for instance [9j) that 
is very helpful in dealing with them combinatorially as well as analytically. Here we will 
represent them interchangeably as two line arrays or as tableaux. A single example of this 
correspondence should be sufficient for our purposes. In the figure below we have on the 
left the two line array, with the list of cars V = (v\,V2, ■ ■ ■ ,v n ) on top and their diagonal 
numbers U = (ui,U2, . . . ,u n ) on the bottom. In the corresponding nx n tableau of lattice 
cells we have shaded the main diagonal (or 0-diagonal) and drawn the supporting Dyck 
path. The component ui gives the number of lattice cells EAST of the i th NORTH step 
and WEST of the main diagonal. The cells adjacent to the NORTH steps of the path are 
filled with the corresponding cars from bottom to top. 
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The resulting tableau uniquely represents a parking function if and only if the cars increase 
up the columns. 
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A necessary and sufficient condition for the vector U to give a Dyck path is that 

(1.2) u\ = and < Uj < + 1 

This given, the column increasing property of the corresponding tableau is assured by the 
requirement that V = (vi,V2, ■ ■ ■ , v n ) is a permutation in S n satisfying 

(1.3) Ui = Ui_i + 1 => Vi > Vi-i 

We should mention that the component Ui may also be viewed as the order of the diagonal 
supporting car Vi. In the example above, car 3 is in the third diagonal, 1 and 8 are in the 
second diagonal, 5, 7 and 6 are in the first diagonal and 2 and 4 are in the main diagonal. 
We have purposely listed the cars by diagonals from right to left starting with the highest 
diagonal. This gives the diagonal word of PF which we will denote o~{PF). It is easily 
seen that a(PF) can also be obtained directly from the 2- line array by successive right to 
left readings of the components of the vector V = (vi, i>2, • • ■ , v n ) according to decreasing 
values of Ui,U2, ■ ■ ■ ,u n . In previous work, each parking function is assigned a weight 

(1.4) W(PF) = t area(PF) g dinv(PF) 

where 

(1.5) area(PF) = u\ + 112 + • • • + u n 
and 

(1.6) dinv(PF) = (x( u i = u j & v i < v j) + x( u i = u j + 1 & v i > v j) 

l<i<j<n 



It is clear from this imagery, that the sum in (1.5) gives the total number of cells between 
the supporting Dyck path and the main diagonal. We also see that two cars in the same 
diagonal with the car on the left smaller than the car on the right will contribute a unit 
to dinv(PF). The same holds true when a car on the left is bigger than a car on the right 
with the latter in the adjacent lower diagonal. Thus in the the present example we have 

area(PF) = 10, dinv(PF) = 4, a(PF) = 31857624, 

yielding 

w{PF) = t l0 q 4 

Here and after, the vectors U and V in the two line representation will be also referred 
to as U(PF) and V{PF). It will also be convenient to denote by VF n the collection of 
parking functions in the n x n lattice square. 

The Shuffle conjecture [10] states that for any partition /j, = (fii, /j,2, ■ ■ ■ , h n we have 
the identity 

(1.7) <Ve„ , ■ ■ ■ h M ) = Yl t are < PF \ din ^ PF h{a{PF) G ^UU^UU • • • UU$) 

PF&VF n 

where V is the Macdonald eigen-operator introduced in [I], e n is the familiar elementary 
symmetric function, h^ 2 ■ ■ ■ hp t is the homogeneous symmetric function basis indexed 
by /i, £1,82, ... ,£t are successive segments of the word 1234 ••• n of respective lengths 
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Hi, fj,2, ■ ■ ■ > fJ>i and the symbol xi a {PF) S fiUJ^LU ■ ■ ■ UU£^) is to indicate that the sum 
is to be carried out over parking functions in VF n whose diagonal word is a shuffle of 
the words £1,82, ■ ■■,££■ In |8j Haglund proved the I = 2 case of (1.7). By a remarkable 
sequence of identities it is shown in [8] that this case is a consequence of the more refined 
identity 
(1.8) 

(A hj E n>k , e n ) = t are < PF \ dm < PF h{a(PF) G £jUJ£ n .j) for 1 < k < n 

PFeTF n +j(k) 

with £j = 12 • • • J, £ n ~j = J + 1 • • ■ J + n, and the sum is over the collection VF n+ j{k) of 
parking functions in the n+Jxn+J lattice square that have k of the cars J + l,...,J + nin 
the main diagonal including car J+n in the cell (1, 1). Here the E n & are certain ubiquitous 
symmetric functions introduced in [3] with sum 

(1-9) E H; i + E n> 2 + • • • + E Ut k = e n 

and Ah- is the linear operator obtained by setting for the modified Macdonald basis in 

im, si- 

(1.10) A hj H^X ]q ,t] = hj[ J2 t^&^H^X-^t] 

(»J)e/i 

More recently, J. Haglund, J. Morse and M. Zabrocki [2] formulated a variety of new 
conjectures yielding surprising refinements of the shuffle conjecture. In |12j they introduce 
a new ingredients in the Theory of parking functions. This is the diagonal composition of 
a Parking function, which we denote by p(PF) and is simply the composition which gives 
the position of the zeros in the vector U = (ui,U2, ■ ■ ■ ,u n ), or equivalently the lengths of 
the segments of the main diagonal between successive hits of its supporting Dyck path. 
One of their conjectures is the identity 

(X7C pi Cp 2 • • • Cp k 1 , /i^j , h^ 2 , • • • , hp Ll ^ = 

(1.11) = t are < PF \ din ^ PF \{a{PF) G £iUU£ 2 UU • • • l±j£t) 

PF<=VF n 
p(PF)=( Pl ,p2,...,p k ) 

being valid for all p \= n and n - Where for each integer a , C a is the operator whose 

action on a symmetric function -F[-X"], in plethystic notation, can be simply expressed in 
the form 

(1.12) C a F[X] = (-l^FlX-^jY^hmlX] 

m>0 



Remarkably, the operators in (1.11) appear to control the shape of the supporting Dyck 



paths. Since in |12j it is shown that we also have the identity 
(1.13) E Ui k = Y^ C Pl C P2 • • • C Pfe 1 

Pl+P2-\ hPfe=n 
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it comes natural to inquire what becomes of Haglund's identity (1.8) when E n j~ is replaced 
by one of the symmetric polynomials C P1 C P2 • • • C Pk 1 . Note, however that since the k 



in (1.8), under the action of A^j controls the number of big cars on the main diagonal, it 
natural to suspect that the combination of A/^ and C Pl C P2 ■ ■ ■ C Pk 1 would result in forcing 
k of the big cars to hit the diagonal according to the composition p = (pi,P2, ■ ■ ■ ,Pk)- 
Miraculous as this might appear to be, computer data beautifully confirms this mechanism 
. . . but up to a point. In fact, following this line of reasoning, one might conjecture the 
identity 

(A/i J C Pl Cp 2 • • • C Pk 1 , e n ) = 

(1.14) Yl t area{PF) q dinv{PF) x(a{PF)el2---JUJj + l---J + n) 

PF€VF J+n (k) 
p(tng(PF))=(pi,p2,...,p k ) 

where p{big{PF)) now refers to the diagonal composition of the big cars, but otherwise the 



sum is over the same parking functions occurring in (1.8). Now that turned out to be false. 
Yet computer data revealed that the following (g-reduced) version of (1.8) is actually true. 
Namely 

<A hj C w C M ---C Pk l, e n )| = t area(PF) x(v(PF) € 12---jUJj+l---j+n) 

p{big{PF))=(p 1 ,p 2 ,...,p k ) 



This circumstance led to the conjecture that (1.14) could be made true by replacing the 
classical parking function "dinv" by a new dinv more focused on the positions of the big 
cars. 

The main result of this paper is a proof of this conjecture. Banking on the intuition 
gained from previous work [7] and using some of the identities developed there with the 
C a and B(, operators we are able to derive the following basic recursion. 

Theorem 1.1. For all compositions p = (pi,P2, ■ ■ ■ ,Pk) we have 

( A/y C pi C p2 • • • Cp fc 1 , e n ^ = t Pl X q k 1 ^A^ j ._ 1 Cp 2 • • • Cp fc B Pl l , e n ^ 
(1-15) + x(pi = 1)(A^C P2 • • • C Pk l , e n -i) 

with B a = wB a w and for any symmetric function F[X] 
(1.16) B a F[X] = F[X - ^]n[zX] 

Now the Haglund-Morse-Zabrocki conjectures also assert that, replacing the C operators 
by the B operators in ( |1.11[ ) has the effect of allowing the controlled Dyck paths to hit 
the diagonal everywhere, including the points forced by the composition p. This led us to 



interpret the first polynomial on the right hand side of (1.8) as a weighted enumeration 
of the collection of parking functions with diagonal word of a shuffle of 12 • • • (J — 1) by 
j ■ ■ ■ n+j — 1 whose big cars hit the main diagonal according to the collection of compositions 
obtained by concatenating (p2, • • • ,Pk) with an arbitrary composition of p\. Guided by this 



A NEW "DINV" ARISING FROM THE TWO PART CASE OF THE SHUFFLE CONJECTURE 5 



interpretation we succeeded to obtain by means of (1.16) a recursive construction of the 
appropriate new dinv and prove the identity 

{^■hjC Pl C P2 ■ ■ ■ C Pk 1 , e n ) = 

(1.17) Yl t are < PF ^q ndin < PF \{a{PF)e\2---JUJj+l---J + n) 

PFeTF J+n (k) 
p(big(PF))=(j>i,p 2 ,...,Pk) 

To carry out all this we need a collection of identities of Macdonald polynomial theory al- 
ready used in previous work. These identities and the corresponding notational conventions 
will be collected in the first section with references to the original sources for their proofs. 
The second section will be dedicated to the proof of Theorem |1.1| All the corresponding 
combinatorial reasoning including the construction of the new dinv is given in the third 
section where our "ndznu" is also given an equivalent somewhat less recursive construction 
with the hope that it may be conducive to the discovery of a direct formula for the new 
dinv which, as in the case of the classical dinv, is closely related to the geometry of the 
corresponding parking function diagram. 
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2. Auxiliary identities from the Theory of Macdonald polynomials 

The space of symmetric polynomials will be denoted A. The subspace of homogeneous 
symmetric polynomials of degree m will be denoted by A =m . We will seldom work with 
symmetric polynomials expressed in terms of variables but rather express them in terms 
of one of the six classical symmetric function bases 

(1) "power 11 {p^}^, (2) "monomiaV {m M } M , (3) "homogeneous 1 '' {h^}^, 
(4) "elementary 11 {e^}^, (5) "forgotten 11 {f^}^ and (6) "Schur 11 {s^}^ . 

We recall that the fundamental involution oj may be defined by setting for the power 
basis indexed by \i = (ni,[j,2, ■ ■ ■ , \~ n 



(2-1) wp„ = {-l) n -% = (-l)W-^ 

where for any vector v = (yi, «2, • • • , v k ) we set \v\ = Y^i=i v i an d K v ) = 

In dealing with symmetric function identities, specially with those arising in the Theory 
of Macdonald Polynomials, we find it convenient and often indispensable to use plethystic 
notation. This device has a straightforward definition which can be verbatim implemented 
in MAPLE or MATHEMATICA for computer experimentation. We simply set for any 
expression E = E(t\,t2, ■ ■ ■) and any power symmetric function 

(2.2) p k [E] = E(t k 1} t k 2 ,...). 

This given, for any symmetric function F we set 

where Qf is the polynomial yielding the expansion of F in terms of the power basis. Note 
that in writing E(t\,t2, ■ ■ ■) we are tacitly assuming that ti,t 2 ,t^, . . . are all the variables 
appearing in E and in writing E(t\, t 2 , . . .) we intend that all the variables appearing in E 
have been raised to their k th power. 



A paradoxical but necessary property of plethystic substitutions is that (2.2) requires 

(2.4) p k [-E] = -p k [E\. 

This notwithstanding, we will still need to carry out ordinary changes of signs. To 
distinguish it from the plethystic minus sign, we will carry out the ordinary sign change 
by prepending our expressions with a superscripted minus sign, or as the case may be, by 
means of a new variables e which outside of the plethystic bracket is simply replaced by 
— 1. For instance, these conventions give for X k = x\ + X2 + • • • + x n 

n 

Pk [--x n ] = (-l)*- 1 ^ 

8=1 

or, equivalently 

n n 

Pk [-eX n ] = -e fc 5>* = (-l)*- 1 ^ 
i=\ i=l 
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In particular we get for X = x\ + X2 + X3 + • • • 

ojpk[X] = p k [-~X]. 
Thus for any symmetric function F £ A and any expression E we have 

(2.5) ujF[E] = F[-~E\ = F[—eE] 
In particular, if F 6 A =k we may also rewrite this as 

(2.6) F[-E] = uF[~E] = {-l) k LoF[E}. 
The formal power series 

(2.7) n = e Xp (Y. P i) 

k>l 

combined with plethysic substitutions will provide a powerful way of dealing with the many 
generating functions occurring in our manipulations. 

Here and after it will be convenient to identify partitions with their (french) Ferrers 
diagram. Given a partition fi and a cell c £ /x, Macdonald introduces four parameters 
I = 1/j,( c )j I' = ^( c ); a = a fi{ c ) an d a! = a^(c) called leg, coleg, arm and coarm which 
give the number of lattice cells of /i strictly NORTH, SOUTH, EAST and WEST of c, (see 
attached figure). Following Macdonald we will set 

1(h) 

(2.8) n(/i) = £Wc) = = E(<-1W 

CG/i CG/^ 2=1 



if coarm ) 



i!5 




(arm ) 



Denoting by \j! the conjugate of /i, the basic ingredients playing a role in the theory of 
Macdonald polynomials are 
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(2.9) n M (g,t)= J] (l-t^V^)), M=(l-t)(l- 9 ) 

ce^t;c^(0,0) 

£>„(?,*) = MB^q,t) - 1 , «J M (g,t) = ]](q a » {c) - t l ^ +1 ){t 1 ^ - q a ^ +1 ), 

together with a deformation of the Hall scalar product, which we call the star scalar 
product, defined by setting for the power basis 

i 

where gives the order of the stabilizer of a permutation with cycle structure \i. 

This given, the modified Macdonald Polynomials we will deal with here are the unique 
symmetric function basis {H^X; q, t)} which is upper triangularly related to the basis 

{sa[j3j]}a and satisfies the orthogonality condition 

(2-10) (H x , = X (A = / uK( (? ,t) 

In this writing we will make intensive use of the operator V defined by setting for all 
partitions fi 

V-ff M = T^H^. 

A closely related family of symmetric function operators is obtained by setting for a sym- 
metric function F[X] 

It is good to keep in mind that, because of the relation e n [Bfj] = for /ihn, the operator 
V itself reduces to A en when acting symmetric polynomials that are homogeneous of degree 
n. 

Recall that for our version of the Macdonald polynomials the Macdonald Reciprocity 
formula states that 

(2 11) H a [l + uD,\ _ H p [l + uD a ] 

(2 - U) rWl-^V) " rWl-ntV) (forallpan-sa,/?) 

We will use here several special evaluations of (2.11). To begin, canceling the common 
factor (1 — u) out of the denominators on both sides of (2.11) then setting u = 1 gives 

H a [MB B ] Hp[MB a ] . 
(2.12) — = — (for all pairs a, p) 

n a Us 



On the other hand replacing u by 1/u and letting u = in (2.11) gives 



(2.13) 



( _ 1)M fi«M = ^mWA (for all pairs a, /?) 
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Since for j3 the empty partition we can take Hp = 1 and Dp = — 1, (2.11) in this case 
for a = \i reduces to 

n-l 

(2.14) H^l-u] = H(l-ut l 'q a ') = (l-tO^-ttfe^-l] 

This identity yields the coefficients of hook Schur functions in the expansion. 

(2.15) H»[X;q,t] = £ s„[X]K x ^q,t) 

Recall that the addition formula for Schur functions gives 

U-u) r {l-u) if^=(n-r,l r ) 
I otherwise 



(2.16) s^l-u] 



Thus (2.15), with X = 1 — u, combined with (2.14) gives for /ihn 

(Hft , S( n _ rjl r)) = e r [B^-l] 
and the identity e r /i n _ r = S( n _ r) ]r) + S( n _ r _ 11 r-i) gives 
(2.17) (H„ , e r h n ^ r ) = e r [B^. 



Since for /3 = (1) we have Hp = 1 and IIg = 1, formula (2.12) reduces to the surprisingly 
simple identity 

(2.18) H a [M] = MB a U a . 

Last but not least we must also recall that we have the Pieri formulas 

(2.19) a) ei H u = J^d^Hf,, b) ei&n = ^ CpLV H u , 

and their corresponding summation formulas (see [2] . [6] . [15] ) 



(2.20) 



(at\(TlT\ k \m h k+l [D^q,t)/tq] if k > 1 
Z^M^W^J " \ B (q )t ) if fe = 



(2.21) EVMKW)' - 

Here i/— >^ simply means that the sum is over v 's obtained from by removing a corner 
cell and //•<— i/ means that the sum is over pi's obtained from f by adding a corner cell. 
It will also be useful to know that these two Pieri coefficients are related by the identity 

(2.22) <V = Mc^ 

Recall that the Hall scalar product in the theory of Symmetric functions may be defined 
by setting, for the power basis 

(2-23) (p x , Pfj.) = Zf, x(A = m), 
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It follows from this that the *-scalar product, is simply related to the Hall scalar product 
by setting for all pairs of symmetric functions /, g 

(2-24) (f,g)^ = (f,^g), 

where it has been customary to let 4> be the operator defined by setting for any symmetric 
function / 

(2.25) 4>f[X\ = f[MX\. 

Note that the inverse of <fi is usually written in the form 
(2-26) f*[X] = f[X/M] 

In particular we also have for all symmetric functions /, g 
(2-27) (f,g) = (f,ug*), . 



The orthogonality relations in (2.10) yield the Cauchy identity for our Macdonald polyno- 
mials in the form 



(2-28) n[-e^] = £ 



which restricted to its homogeneous component of degree n in X and Y reduces to 



(2.29) = •£ "" {X y Y] 



Note that the orthogonality relations in (2.10) yield us the following Macdonald poly- 
nomial expansions 
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Proposition 2.1. For all n > 1 we have 



fjhn 



W, 



(2.30) 



fi\-n 

e) ei[X/M] ? 



fi\-n 



/) 



E 



Finally it is good to keep in mind, for future use, that we have for all partitions [i 
(2.31) T^H^X-l/q^/t] = H^[X;q,t\ 

Remark 2.2. It was conjectured in [5] and proved in [11] that the bigraded Frobenius 
characteristic of the diagonal Harmonics of S n is given by the symmetric function 

T„HJX; q, t)MBJq, t)UJq, t) 



(2.32) 



DH n [X;q,t] = £ ' " 

fi\—n 



Surprisingly the intricate rational function on the right hand side is none other than V6 n . 
To see this we simply combine the relation in (2.18 ) with the degree n restricted Macdonald- 



Cauchy formula (2.29|) obtaining 
(2.33) 



e n [X] = e n [™ 



M 



E 

fi\-n 



Wi, 



This is perhaps the simplest way to prove (2.30) f). This discovery is precisely what led to 
the introduction of V in the first place. 
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3. Proof of the basic recursion 



To establish Theorem 1.1 we need some preliminary observations. To begin we have the 
following reduction. 

Theorem 3.1. For all p = (p±,P2, ■ ■ ■ ,Pk) |= n & n d j > we have 

{^hjC pi C P2 ■ ■ ■ C Pk l , e n ^> = t Pl 1 ^A^ j ._ 1 Cp 2 • • • Cp fc B Pl l , e n ^ 
(3-1) + x(pi = 1)(A^C P2 • • • C Pk l , e n _i> 

if and only if , with C* and B* the *- scalar product duals of C a and B a , we have 
(3.2) C* a A hj h n [§] = t a - 1 B* a A hj _ 1 h n [^} + x(a = l)A hj h n ^ 1 [§] 

for all j > and 1 < a < n. 

Proof. It is shown in [12J that the operators C a and B& satisfy the commutativity relations 

qC a B b = B b C a (for all a, b > 1) 



Using these identities, (3.1) becomes 



{^hjC Pl C P2 ■ ■ ■ C Pk l , e n ^ — t pi ^A^_ 1 Bp 1 Cp 2 • • • Cp fc l , e n ^ 

+ x(Pi = l)(A/ l ,C P2 • • • C p . 1 , e. 



--ra— 1 , 



Passing to *-scalar products, and using the identity in (2.27), we can next rewrite (3.1) in 

the form 

(3.3) 

(A h .C Pl F[X] , h* n )^ = ^(A^Bp.FlX] , h* r X + X (P1 = l){A hj F[X] , 

and the validity of this identity for every symmetric function F[X] that is homogeneous 
of degree n — p\ is equivalent to (3.1 ) since when P2, ■ ■ - Pk are the parts of a partition the 



polynomials C P2 • • • C Pk l are essentially elements of the Hall-Littlewood basis. Now, since 



all the operators Ap are self adjoint with respect to the *-scalar product, (3.3) in turn can 
be rewritten in the form 

(F[X] , C* n A h .h* n )^ = t^-\F[X] , B^A^K)^ + X (P1 = 1)(F[X] , A^i), 



and this identity (for all p\ > 1) is equivalent to (3.2) due to the arbitrariness of ^[-X]. 
This completes our proof. □ 



Our next goal is to prove (3.2). To begin we have the following auxiliary identity. 
Proposition 3.2. 



(3.4) 



A hjh n [w] = E^^](-1) S E^4^(M-^)]" 



s=0 



v\-s 
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Proof. Using (2.30) c) and the definition of the operator A/,, we get 



^m§] = y.'"":: x >>r : 



w. 



^ w, J 1 



M + mJ 



Ill l-MiWl 



M 



j 



if. m 

E. rj_] v - ^ W T U H V [MB^ — 1] 

s=0 fihn v\-s 

4^h riiv T 'v WW HAMB^ - i] 



s=0 

j 



(by the reciprocity in (Q) = E S(-l)^^ ^3^g^ ~ 1] 



s=0 
3 



W, 



T, 2 



(ty@)) = EW£]£:f(-irV*(*"-£)] 



s=Q 



This gives Q since e n [X(5, - ^)] = (-l)"/i n [X(^ - £„)] . 



It is good to keep in mind that in particular we also have 



(3.5) 



s=0 



(3-6) & hj -M§\ = D-^-i-.^E^Wi--^)] 



s=0 



Next we have 
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Proposition 3.3. 

(3.7) 

C* a A hj h n [§) - x (a = l)A hj h n - 1 [ 

3-1 



M \ 



T. 



+ e- 1 £ h^si^K-iy E 77 E h — ifS h — ~ ^)] h - i 1 ~ MB A 



s=0 



t\-s v=a 



Proof. The identity in (|3.4|) gives 
(3.8) 



C* a A h .h n [ x 



T 



All 



E h i-siU-iy E t c « hn ^ ~ B ^ ■ 



uhs 



Now it was shown in [7] that for all P[X] £ Awe have 

c* a p[x] = (=i r -ip[x-*-f]n[^ x 

This gives 



*)- 



(- q y-i C * a h n [X(±-B v )] = h„[(X Bv)]n[ 



-tzX - 
9(1-*)- 



q(l-t). 



r=0 



£ h n _ r [x(i - b u )] {-\) r h r [mb u - 1] n [ 



-ezX - 
q(l-t). 



r=0 



cX 



£ K-r [X{ij ~ B v )) {-Ifhr [MB U - 1] h r „ a 

r=0 
n 

hn-r [X(jj - B v )} { -^h r [MB V - 1] hr-a [j 



■$] 



or better 

C*/t n [X( M - Bj,)] 



^ /i n _ r [X(jj - B v )] -kxh r [MB U - 1] h r _ a [= 



and the last sum in (3.8) becomes 
(3.9) 



E «r^-a[^] E :f V-r[*(£ - ^)] /i r [mb„ - 1] 

Mil 
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Using the summation identity in (2.20) in the form 

h r [MB U - 1] = (tqy- 1 ]T Mcw^r 1 - X (r = 1) 

T — 

we get 



+ E ^ [fS\ E „f ^ - B v )] X (r = 1) 



v\~s 



T 2 



E^-H] E - ^)] E Mcur{%r x 



v\~s 



T, 



+ X {a=l)Y J — h n-i[X(i I -B v 



u\-s 



r+l 



rhs-1 i>«-t 



T 2 



+ x (a = l)^^ n _ 1 [X(^-^ 

. vhs 



Using (2.22) and the fact that there are no partitions of size —1, ( |3.8[ ) becomes 

3 n rp2 

c:A hjhn [§] = -^Ev^K-^E^t^f] EfE^^F-^lMir 



+i 



8=1 



W r 

rhs-1 fc"*-T 



T 2 



+ X (a = i) E hj-siM- 1 ) 3 E ^r^-i - B »j\ 



s=0 



and (3.5) gives 



C* a A h .h n [ 
(3.10) 



^1 

Mi 



Tl 



s=l 



rhs-1 I>«-T 



-T 

+ x(a = l)A ftj /i„_i[^] 



r+l 



Since B v = B T + ip^ we derive that 

n— r 

Y,K-r[X{jg ~ B v )}d„ T {^y +l = Y, h n-r-u[X{ij-B T )]h u [-X]Y,dur{^r 

T U=0 Vr- T 

n 

= Y, h n-v[x(^-B T )]h v „ r [-x]Y,d„A%: 



.r+u+1 



v+1 



Hi 



A. DUANE, A. M. GARSIA AND M. ZABROCKI 



Using the summation formula in (2.21) in the form 



h k -i[l-MB u ] ifjfe>l 
1 ifk = . 



together with the fact that v > r and in ( |3. 10 > we have r > a > 1 we obtain 

= h v [l-MB T ] 



and (3.10) becomes 
C* a A hj h n [§] - x(fl = l)A hj h n - 1 [ 



x_- 

M. 
n 



8=1 



rha-1 



\u> r 



E ^ - B t)] hv-r [ ~ X] h v [1 - MB T 



v=r 
3 



-^E^iiK-D* E ^EEM^ 



8=1 



rhs— 1 ti=a r=a 



/l r _ a [=^] /in-, - B T )] fc„ [1 - MB T 



^E'y-.tiK-i)* E -f E^-^K-^m-^M 1 -^ 



8=1 



10. 

rhs— 1 f=a 



This best rewritten in the form 
CiA/y/^ff] - X {a = l)A hj h n - l [§] = 

= t- 1 e (-ir E ^ E ^ - ^)] ^ t 1 - MB r] 



s=0 



rhs v=a 



proving (3.7) and completing our proof of Proposition 3.3 



□ 



Let us now work on B*A/ l ._ 1 /i n [jg>]. Here we use the identity in (3.6), that is 



T 2 



Mi 



EVi- s [^](-i) 8 E^^[^(f-^ 



s=0 



and the identity 



b:p[x] = P [x + M.]n[=g- 
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(established in [7]) that gives the action of the operators B* to obtain 



E ;f b :M*(s - ^)] 



s=0 

i-i 



r, 



E(-!) s Vi- [£] E :f M(* + t )(m - ^)] n [t 



zX] 



T 2 



s=0 uhs 

EC-^-i- [»] E E ^ " ^ [ MB * - !] n ft 

s=0 

j-1 



vhs r=0 



T, 



Ei- 1 ) 8 ^ [if] E E ^ - ^ E 1 - M ^ l*S 



s=0 

and this may be rewritten as 
i-i 



B*A 



r, 



£>ir vi- 8 [&] E E lxffl ^ t 1 -^ 



s=0 



Comparing with the right hand side of ( |3.7[ ) we see that we have established the identity 

Ml 

This completes our proof of (3.2) an consequently also the proof of Theorem 



e- 1 BtA h ._ 1 h n [§] = C* a A hj h n [§] -x(a = l)A h .h n - l [ 
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4. The construction of the new dinv 



Let VF(J, n) denote the collection of Parking functions on the J+n x J+n lattice square 
whose diagonal word is a shuffle of the two words £ j = 12 . . . J and £j^ n = J + 1 . . . J + n 
with car J + n in the (1,1) lattice square. In symbols 



Before we can proceed with our construction of the new dinv, we need some preliminary 
observations about this family of parking functions. To begin we should note that the 
condition that the diagonal word be a shuffle of 12 • • • J with J + 1 • • ■ J + n, together with 
the column increasing property of parking functions, forces the columns of the Dyck path[j] 
Where by the length of column i of a Dyck path D we refer to the number of NORTH steps 
of D of abscissa i supporting a PF G VF(J,n) to be of length 2 at most. The reason for 
this is simple: as we read the cars of PF to obtain a(PF) from right to left by diagonals 
starting from the highest and ending with the lowest the big cars ( J + 1, . . . , J + n) as well 
as the small cars (1,2,..., J) will be increasing. Thus we will never see a big car on top of 
a big car nor a small car on top of a small car. So the only possibility is a big car on top 
of a small car, i. e. columns of length 2 at most as we asserted. 

This yields an algorithm for constructing all the elements of the family VF(J, n). Let us 
denote by "red(PF)" , and call it the "reduced tableau " of PF, the configuration obtained 
by replacing in a PF G VF(J,n) all big cars by a 2 an all small cars by a 1. We can 
simply obtain all the reduced tableaux of elements of VF(J, n) by constructing first the 
family T>j^ n of Dyck paths of length n + J with no more than J columns of length 2 and all 
remaining columns of length 1. Then for each Dyck path D 6 T>j^ n fill the cells adjacent 
to the NORTH steps of each column of length 2 by a 1 under a 2, then fill the columns of 
length 1 by a 1 or a 2 for a total of J ones and n twos. 

Clearly each PF G VF(J, n) can be uniquely reconstructed from its reduced tableau by 
replacing all the ones by 1, 2, . . . , J and all the twos by J + 1 ... J + n by diagonals from 
right to left starting from the highest and ending with the lowest. It will also be clear 
that we need only work with reduced tableaux to construct our new dinv. However, being 
able to refer to the original cars will turn out to be more convenient in some of our proofs. 
For this reason we will work with a PF or its red(PF) intercheangeably depending on the 
context. 

This given, we have the following basic fact 
Proposition 4.1. For any 



(4.1) 



VF(J, n) 



{PF e VF J+n : a{PF) £ £jUJ£ n -j & J + n G (1, 1)} 



PF 



G VF(J, n) 



Ul U 2 



if we set {i% < 12 < ■ ■ ■ < ik} 



{i G [1, J + n] : > J} then the vector 



Ub(PF) 



X*) 
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gives the area sequence of a Dyck path, which here and after will be referred to as the Dyck 
path "supporting" the big cars of PF. 

Proof. Since car J + n is in the (1,1) lattice square it follows that =0. Thus we need 
only show that 

u i s < u i s -i + 1 f° r all 2 < s < /c 
By definition vi s _ l and Vi s are two successive big cars in PF that means that for i s -\ < j < 



the car Vj is small and thus, except perhaps for j 



1 , the car Vj must be in a column 



of length 1. In particular we see that we must have Uj < Ui s _ 1 for, the first violation of this 
inequality would put a small car above a big car (for j = i s _i + 1) or a small car above a 
small car for j > i s _\ + 1. This gives Ui s < Ui g -i + 1 as desired with equality only if car 
Vi s is at the top of a column of length 2 and all the small cars in between Vi s _ 1 and Vi s are 
in the same diagonal as Vi g _ 1 . □ 

In view of this result we are now going to focus on the subfamilies 7~*F j(p) of *PF 
consisting of its elements whose big cars have a supporting Dyck path which hits the 
diagonal according to a given composition p = (pi,p2, • • • ,Pk) 1= n - O ur goal is to construct 
a statistic "ndinv " which yields the equality. 

(4.2) (A hj C pi C pi ■ ■ ■ C pk l , e n ) = £ t area(PF) q ndinv(PF) 

PFeVFj(p) 

But before we do this it may be good to experiment a little by constructing some of these 
families. 



By reversing the argument we used to prove Proposition 4.1 we can start by constructing 
all the Dyck paths with the given diagonal composition then ad all the cars as required 
by the definition of a family. This is best illustrated by examples. Say we start with with 
p = (3, 2). In this case there are only two possible Dyck paths as given below on the left . 



















to 




(2) 










n 










(2) 



























® 






2! 


















(2) 











On the right we added the 2's and their corresponding diagonal numbers. Now the least 
number of l's we need to add to get a legal reduced diagram is 3 for the first and 2 for the 
second as shown below 



(4.3) 
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Now a MAPLE computation yields the polynomials 

(4.4) (A fca C 3 C 2 l , e 5 ) = t 3 q 4 
and 

(4.5) <A fe3 C 3 C 2 l , e 5 ) = t 3 (q 4 + q 5 + q 6 + q 7 + q S ) + t\q A + q 5 + q 6 ) + t 5 q 4 



To compute the classical weight t area q dinv it is better to have a look at the non-reduced 
versions of the two tableau above. Namely 



(4.6) 



















































































































00 

















and 

















































(4; 












(5 


























(J] 















Now in the first PF, the pairs (3,7), (1,5) and (6,2), are the only ones contributing to 
the dinv and the sum of the area numbers is 5, so its classical weight is t 5 q s . Similarly, 
the pairs contributing to the dinv on the PF on the right are (2, 6), (5, 1) and (4, 1) and 
the area numbers add to 3, so its classical weight is t s q 3 . The latter is not the same as 
what comes out of (4.3). The area is OK but the dinv in not. The calculation in (4.4) 



thus asserts that the "new dinv" should be 4. Similarly, as we will show in a moment, the 



calculation in (4.5) yields that the new dinv of the PF on the left of (4.6) should be 4 
again. In fact, it turns out that none of the 8 parking functions we obtain by inserting an 



extra 1 in the reduced tableau on the right of (4.3) have area 5 thus the last term in (4.5) 
can only be produced by the PF on the left of (4.6). We give below the 8 above mentioned 
reduced tableaux with the extra 1 shaded 



(4.7) 



44 















'2) 
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Therefore the reduced tableaux of the family VF '3 ([3, 2]) are 9 altogether, as predicted by 
(4.5), namely the 8 above together with the tableau on the left of (4.3). Computing their 
classical weight and summing gives 

t area(PF) q dvnv(PF) = ^4 + ^5 + ^6) + ^ (g 3 + g A + g5) + ^3 

PFeVF s ([3,2]) 



the first 8 terms from (4.7) and the last from the left tableau in (4.3). As we see this is 
not quite the same polynomial as in (4.5). Note that the area again works but the classical 
dinv does not!. 

For a while in our investigation this appeared to be a challenging puzzle. The discovery 
of the recursion of Theorem 1.1 completely solved this puzzle but, as we shall see, it created 
another puzzle. 



Let us have a closer look at (1.16), namely the identity 



(A hj C Pl C 



V2 



C Pk 1 , e n/ 



v /i,/_iCp 2 



Cp^Bpj 1 , e n ) 



(4.8) 



+ x(Pi = 1 )( A /i i C P2 ■ ■ • C Pk l , en^) 



Setting for a composition p = (pi,P2, ■ ■ -Pk 



n 



j.area(PF) qndinv(PF) 
PFeVFj(p) 



our conjecture, together with the identity 



B P1 1 



-pi 



(<?i><72,---,<?f)hpi 



C<?1 Cq 2 ■ ■ ■ C Ql 1 



proved in [TU], translates (4.8) into the recursion 

(4.9) Iij([ Pl ,p 2l ...p k ]) = t^- 1 q k - 1 Y,^J-i([P2,---Pk,r]) + X (pi = mA\P2,...Pk}) 

r\=pi 

where the symbol \p2, ■ ■ ■Pk, v] represents the concatenation of the compositions (p 2 , ■ ■ -p k ) 
and r. This strongly suggests what should recursively happen to the new weight of our 
parking functions by the removal of a single (appropriate) car. That is, if the chosen car is 
small there should be a loss of area of p\ — 1 and a loss of ndinv of k — 1 and if the chosen 
car is big no loss of any kind. 

Starting from this observation and further closer analysis of (4.8) led us to the following 
recursive algorithm for constructing "ndinv" . 

This is best described by working with the corresponding reduced tableaux. To begin 
it will be convenient to start by decomposing each red(PF) into sections corresponding 
to the parts of the given composition. To be more precise it is best viewing our two line 
arrays as unions of vertical dominos. For instance the red(PF) below, which is none other 
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than the minimal ones obtained from the Dyck path on the right 



2 o 
1 . 



. 2 o o 

, 1 o , 

2 a . . 

1 . . . 



, 2 o o 

. 1 o . 

2 o . . 

1 . . . 



will be viewed as the sequence of dominos 
(4.10) 





2 




1 




2 




1 




2 




2 




1 




2 




1 




2 




2 




1 




2 
















1 




1 




2 














1 




1 




2 














1 





Thus the corresponding PF belongs to the family VF§{$, 3, 2]) and as such will be divided 
into 3 sections, one for each part of [3, 3, 2]. To do this we simply cut the sequence in (4.10 ) 
before each domino [?] obtaining the three sections 



(4.11) 

Since here p\ 





2 




1 




2 




1 




2 








2 




1 




2 




1 




2 








2 




1 




2 
















1 




1 




2 


















1 




1 




2 


















1 





3 > 1, 



(|4.8|) suggests that we should remove a 1 from the first section, 

k 



then process it somewhat to cause a loss of dinv of 2 ( = k — 1), and loss of area 2 
(= pi — 1). Taking a clue from the classical dinv, we can see that the first small car in 
the corresponding PF would contribute a unit to the classical dinv with the big cars to 
its right in the main diagonal. The latter of course correspond to the dominoes [q] that 
begin each of the following sections. Thus the desired loss of dinv can be simply obtained 



by bodily moving the first section to the end, and removing the [J] obtaining 

2 




(4.12) 



We may thus view that the removed domino contributed a unit to dinv for each domino 
to its right. But we still have not accounted for the loss of area and worse yet we will 



L0- 

now have a big car on top of a big car. Since (4.8) tells that the loss of area should be 
pi — 1 then it must be equal to the number of big cars in the moved section, minus one. 
This means that we can fix both problems by making the domino replacements [J — > [q] 



and 



], obtaining 



(4.13) 
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but that creates a new problem, since the succession [J would put a small car on top 
of a big car. We will fix this final problem by simply switching the 1 with the 2 obtaining 



(4.14) 



which gives the domino sequence of the red(PF) below 



. 2 o 

2 o . 
1 . . 



2 o a 
1 o . 



on the right of which we have depicted the Dyck path supporting the big cars. 

In the case that p\ = 1 there will be only one big car in the first section and if there 
are small cars they all must be on the main diagonal. In this case we can process the first 
section as we did for pi > 1. If there are no small cars then the first section consists of the 



single domino L] and (4.8) suggests that we should simply remove it with no further ado. 



To carry out our definition of "ndinv" rigorously and in full generality, we will break our 
argument into three separate steps. In the firs step we use the ideas stemming from the 
above example to construct a bijection 
(4.15) 

|J VFj^ 1 ([p 2l ... Pk ,q\) if Pi > 1 

q\=pi 

VF J ^ 1 ([p 2 ,...p k ,l}) VFj(\p 2 ,...p k ]) if pi = 1 

In the second step we define "ndinv" by setting for each PF E VFj(\pi, P 2, ■ ■ -Pk\) 

' k - 1 + ndinv ($(PF)) if J > 
if J = 



$ : VFj([p 1 ,p 2 ,...pk]) 



(4.16) 



ndinv(PF) 



From step 1 and step 2 it will follow that the polynomials Hj([pi,p2, ■ ■ ■ ,Pk\) satisfy the 
same recursion as the polynomials (A/ l7 C pi C P2 • • • C Pk l , e„'). 



In the third step we establish the equality in (4.2) by verifying the equality in the base 
cases. 



In our first step, starting with a PF E PFj(\pi,p2, 
following procedure. 



■ Pk\) we construct $(PF) by the 



Cut the domino sequence of red(PF) into sections starting at the dominos 
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(1) If the first section does not contain a domino [J] 

• remove its only domino [jj] from the sequence of dominos 

(2) If the first section contains a domino [q], work on the first section as follows 

• Remove its first domino [J] 



2 a ] make the replacement [ ] — > [ Q 2 ] 
2 nil i j ■ , i r 2 nil . r 1 u2i 



• For each (but the first) domino [ ] make the replacement L , [(( _| 



• i/ adjacent pairs [ a _J [ a ] are created make the replacements [ a _ li[ai , L a _uL a j 

• Cycle the modified first section to the end of the sequence of dominos 

In any case we let PF' be the parking function corresponding to the resulting domino 
sequence and set 

§(PF) = PF'. 



It is clear that <1> maps the left hand side of (4.15 ) into the right hand side. To show that 
$ is a bijection we need only show that the procedure above can be reversed to reconstruct 
PF from PF' for any PF' in the right hand side of (4.15). We will outline the salient 
steps of the reversed procedure. 

Note first that since our target PF = is to be in VFj(\pi,p2, ■ ■ -Pk]) we 

already know the diagonal composition of the Dyck path of the big cars of PF. Thus we 
can proceed as follows 

(1) Say PF' £ VF j([p2, ■ ■ -Pk]) (that will only occur when pi = 1) 

• Then PF is the parking function obtained by prepending f 2 ] to the domino 
sequence of PF'. 

(2) Say PF' G VF j_i([p2> • • -Pk-, 1]) (that will only occur when p\ = 1) 

• Then PF is the parking function obtained by inserting [J] immediately after 
the first [q] in the last section of PF' , then cycle back the last section to be the 
first in the domino sequence. 

(3) Say PF' G VFj^{\p 2 , . . . Pk , q]) for a q \= Pl - 1 > 

• Let last(PF') be the domino sequence obtained by removing the first k — 1 
sections from the domino sequence of PF' . 

• Modify last(PF') by inserting a [q] immediately after its first [jj]. 

• For a > 1 replace, in last(PF'), each pair [^^[a] by the pair [ 2 ][ ]. 

(note that for this to put a big car on top of a big car we must have a [ a \] 
preceding the [^J, but that [ (l 2 1 ] will also be replaced either by this step or 
by the next steps) 

• For a > 1 replace each [ 2 ] preceded by a [ a ] in last(PF') by [ a ^_ 1 ] 

• Replace each [ 2 ], except the first by a [ 2 ] 

(note if a replaced [jj] is preceded by a [q] then that [ ] itself will also be 
replaced by [ 2 ]) 

• The modified last(PF') followed by the the first k — 1 sections of PF' gives 
then the domino sequence of our target PF. 
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This completes our proof that <3? is bijective. 

Since moves EAST, by one cell, p\ — 1 big cars it causes a loss of area equal to p\ — 1. 



Thus the definition in (4.16) combined by the bijectivity of proves the recursion in (4.9). 

It remains to show equality in the base cases which, in view of the definition in (4.16) 
should be characterized by the absence of small cars. 

Now it easily seen, combinatorially, that VFo(\p]) is an empty family except when all the 
components of p are equal to 1. To see this note that it is only the presence of small cars 
that allows the supporting Dyck path of one of our PF's to have columns of lengths 2. But 
if all the columns are of length 1, the "area" statistic is and the Dyck path supporting 
the big cars can only have area sequence a string of O's. But in this case the family reduces 
to a single parking function which consists of cars 1, 2, • • • , n placed on the main diagonal 



from top to bottom. Thus it follows from our definition of Hj(p) and (4.16) that 

I if some pi > 1 



IIo([pi,P2, • • • ,Pk}) 



1 if all pi = 1 



Since by definition A/j reduces to the identity operator, the equality for the bases cases 
results from the following fact 

Theorem 4.2. For p = (pi,P2, ■ ■ ■ ,Pk) |= n we have 

> I if some Pi > 1 



(4.17) <C Pl CV--C Pfe l 



if all pi = 1 



Proof. Recall from (1.2) that for any symmetric function F[X] we have 



(4.18) C a F[X] = (-lr- 1 F[X- 1 ^]^2z m h m [X] 

m>0 



In particular it follows that for any Schur function s\ we have 

C a s x [X] = (-ir-^SA/^^Vs-lK+MW 

This gives for a + |A| = n 

(4.19) <C aSA [X] , e n ) = {-IT' 1 Y, s„[l/q ~ l^A/A+H > e ™> 

and the Littlewood-Richardson rule gives 

(sx/^K+M , en) = (s x /im , K+\n\ en ) = 



unless a + \fi\ = 1, Thus for a > 1 (4.19) reduces to 



r^ n \ , ,,. i, , > ll if a = 1 and A 

(4.20) <<wx],«.) = (-ir'<.*»..«.> - | otherwise 

Since 

c a i = Hr 1 ^ 
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the first case of (4.17) follows immediately from (4.20). On the other hand even when all 



the pi are equal to 1 in successive applications of Ci only the term corresponding to si* 
in the Schur function expansion of C™1 will survive in the scalar product 



(CiC™l , e m+ i) 



since from (4.20) it follows that Cfl 



1, the second case of (4.17) is also another 



consequence of (4.20). 



in (1.17). 



This completes of proof of (4.17). This was the last fact we need to establish the equality 

□ 



Remark 4.3. As we already mentioned, our definition of ndinv creates another puzzle. 
Indeed, the classical dinv can be immediately computed from the geometry of the parking 



function or directly from (1.6) which expresses it explicitly in terms of the two line array 



representation. For this reason we made a particular effort to obtain a non recursive 
construction of ndinv and in the best scenario derive form it an explicit formula similar 



to (1.6). However our efforts yielded only a partially non-recursive construction. In our 



original plan of writing we decided to include this further result even though in the end 
it yields a more complex algorithm for computing ndinv than from the original recursion. 
This was in the hope that our final construction may be conducive to the discovery of an 
explicit formula. It develops that during the preparation of this manuscript a new and 
better reason emerged for the inclusion of our final construction. It turns out that Angela 
Hicks and Yeonkyung Kim have very recently succeeded in discovering the desired explicit 
formula by a careful analysis of the combinatorial identities we are about to present. The 
results of Hicks-Kim will appear in a separate publication [T3] , 



For our less recursive construction of ndinv it will be convenient to make a few changes 
in the domino sequences. To begin, we shall use the actual car numbers at the top of the 
dominos rather than 1 or 2. We do this, so that we may refer to individual dominos by 
their car as the corresponding area number on the bottom is being changed. But now, to 
distinguish big cars from small cars we must in each case specify the number J of small 
cars. Secondly, we will have sections end with a big car, rather than begin with a big 
car. This only requires, moving the initial big car to the end of the domino sequence. For 
example, the parking function below whose domino sequence was given in (4.10) has J = 5 
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thus cars 1, 2, 3, 4, 5 are small and 5, 6 . . . , 13 are big. 



(4.21) 
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its domino sequence is now 
(4.22) 
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and its decomposition into sections is as shown below 
(4.23) 



[2 




For convenience we will now need to use the symbols "[*]" and "[*]" to respectively 
represent "big" and "small " car dominos. To motivate our second construction of ndinv 
we will begin by modifying our first construction to adapt to these new domino sequences. 

(1) The recursive construction will now consist of as many steps as there are dominos 
in the domino sequence 

(2) At each step the first domino of the first section is removed 

{a) when we remove an [ S Q \, the section is cycled to the end after it is processed as before 
(b)when we remove a [jj], it is because the section consisted of a single big car domino. 

(3) The removal of an [*] contributes to ndinv the number of[^]'s to its right minus one. 
There are a few observations to be made about the effect of the cycling process. To 

begin note that when the domino sequence consists of a single section, no visible cycling 
occurs. However, even in this case, for accounting purposes, it is convenient to consider 
all of its dominos to have been cycled. With this provision, each domino in the original 
domino sequence will have an associated cycling number c that counts the number of times 
it has been cycled before it is removed. 

Based on these observations, a step by step study of our recursive construction of ndinv 
led us to the following somewhat less recursive algorithm. It consists of two stages. In the 
first stage, the domino sequence is doctored and wrapped around a circle to be used in 
the second stage. The second stage uses circular motion to mimic the cycling of sections 
that takes place in the recursive procedure. To facilitate the understanding of the resulting 
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algorithm we will illustrate each stage by applying it to the parking function in (4.21). 
More precisely we work as follows 

Stage I 

• Move each [*] in the domino sequence a places to its left and increase the area 
number by 1 of each domino [ b ] that is being by-passed 



(for instance the domino sections in (4.23) become) 



(4.24) 



12 




Next wrap the resulting sequence clockwise around a circle with positions marked 
by a "o " 



(4.25) 





(we also place a bar "|" to separating beginning and ending dominos ) 
(the o's will be successively changed to »'s during the second stage) 

Stage II 

• set ndinv = and set the auxiliary parameter c to 1 

• mark the first domino by changing its "o " to a "• ". 

• cycling clockwise from the first domino to the bar find the first call it "endsec" 

• cycling clockwise from endsec to the bar add 1 to ndinv each time we meet a [?]. 



on the right in (4.25) we have darkly boxed the first domino and the endsec) 



(and lightly boxed the two ndinv contributing big car dominos) 
While there is a domino that has not been marked repeat the following steps 

• cycling clockwise from the last endsec mark the first unmarked domino 

• If in so doing the bar is crossed add 1 to c 

If the domino is a [ s ] then clockwise from it find the first [ b ] with a < c, call it "endsec" 
then cycle clockwise from endsec back to this [ s ] 
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• for each encountered unmarked [^] add 1 to ndinv provided a < c if the bar is not 
crossed or a < c + 1 after the bar is crossed 

(the desired value of ndinv is reached after all the the small car dominos are marked) 




(4.27) 




(thus, in this case, ndinv = 14, which is the total number 
of lightly boxed dominos in the previous five configurations ) 
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Remark 4.4. We will not include a proof of the validity of this second algorithm, since A. 
S. Hicks and Y. Kim, using their discoveries, are able to provide in [13] a much simpler and 
more revealing validity argument than we can offer with our present tools. Here it should 
be sufficient to acknowledge that the auxiliary domino sequence resulting from Phase I 
together with the c statistic constructed in Phase II have ultimately been put to such 
beautiful use in subsequent work. 

Before closing we should note that our ndinv may have an extension that can be used 
in a more general settings than the present one. To see this, let us recall that the 2 part 
case of the Shuffle Conjecture, proved by J. Haglund in jS], may be stated as follows 

(Ve n , hjhn-j) = 

(4.28) Yj t are < PF ^q din ^ PF \{a{PF) 6 12 • • • jUJj + 1 • • • n) 

PFeVF n 



Now replacing n by n + 1 — j and J by j in (1.17), for (pi,P2, ■ ■ ■ ,Pk) \= n + 1 — j we get 
(AhjC Pl C P2 • • • C Pk 1 , e n+ i_j) = 

(4.29) t area{ - PF) q ndinv{PF) x{a{PF) € 12 • • • jUJj + 1 • • • n + 1) 

PFeVF n+1 (k) 
p(big(PF))=(pi,p 2 ,...,Pk) 

This given, since is was shown in |12j that we may write 

Cn+l = C pi C p2 • • • Cp l(p) 1 

p\=n+l 

it follows, by summing ( 4.29| ) over all compositions of n + 1, that we also have 
(4.30) 

(A hj e n+1 ^ , e n+1 ^) = ]T« t area(PF) q ndinv(PF) x{a ( pF) £ i2 • • • jUUj + 1 • • • n + 1) 

PFdVF n+1 

where the "(*)" is to signify that the sum is over all parking functions in the n + lxn+1 
lattice square which have the biggest car n + 1 in cell (1, 1). But it was also shown in [8] 
that we have 

\^hj e n+l-j > e n+ i_j) = (Ve„ , hjh n _j) 

Thus (4.30) may also be rewritten in the form 
(4.31) 

<Ve n , hjhn-j) = Y, ( * ] t area(PF) g Wira(PF) x(^(^) G 12 • • • jUUj + 1 • • • n + 1) 

PFdVF n+1 

which gives another parking function interpretation to this remarkable polynomial. It is 
natural then to ask if this kind of result involving the same ndinv, or a suitable extension 
of it, may give a new parking function interpretation to any of the polynomials occurring 
on the left hand side of (1.7). If that were the case then that would provide an alternate 
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form of the Shuffle Conjecture. It is interesting to note that computer exploration has led 
us to conjecture that for p = (pi,P2, ■ ■ ■ ,Pk) \= n the polynomials 

{^■hj 1 ej 2 C pi C pi • • • C Pl l , e n ) 

have non negative integer coefficients. This yields us yet another avenue by which the 
results of this paper can be extended. It should be worthwhile to pursue these avenues 
in further investigations on the connections between Parking Functions and the Theory of 
Macdonald Polynomials. 
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